Time-dependent localized Hartree-Fock potential 
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By minimizing the difference between the left- and the right-hand sides of the many-body 
time-dependent Schrodinger equation with the Slater-determinant wave-function, we derive a non- 
adiabatic and self-interaction free time-dependent single-particle effective potential which is the 
generalization to the time-dependent case of the localized Hartree-Fock potential. The new poten- 
tial can be efficiently used within the framework of the time-dependent density-functional theory 
as we demonstrate by the evaluation of the wave-vector and frequency dependent exchange kernel 
f x {q,oj) and the exchange shear modulus fi x of the homogeneous electron gas. 



Time-dependent (TD) density-functional theory 
(DFT) fl is in the perpetual search of effective single- 
particle potentials which accurately (in ideal - exactly) 
map the propagation of an interacting many-body 
system onto that of the non-interacting one. While 
the exchange-correlation (xc) functionals based on 
the Local-Density Approximation (LDA) 0, Q and 
its semi-local refinement of the Generalized Gradient 
Approximation (GGA) 4-6] have proven very successful 
in the ground-state DFT, in TDDFT the usefulness of 
(semi-) local approaches is limited due to the fundamen- 
tal spatial non-locality of the exact time-dependent xc 
functional 0- 

Beyond LDA and GGA, the concept of the opti- 
mized effective potential (OEP) [1, 0] plays one of the 
key roles in the systematic non-heuristic construction 
of DFT. In the ground-state case, OEP is defined as a 
single-particle potential which minimizes the many-body 
Hamiltonian expectation value on the Slater-determinant 
wave-function. From the DFT perspective, OEP is the 
first term in the adiabatic connection series in the powers 
of the interaction constant (exact exchange) [loj ]. The 
generalization of the OEP to the time-dependent case 
has been proposed which, however, is impractical 
for applications due to the formidable complexity of the 
integral equation involved. 

In this work we propose and implement an alterna- 
tive approach to the development of the time-dependent 
single-particle effective potential for many-body prob- 
lems. This is based on the variational principle of the 
minimization of the difference between the left- and right- 
hand sides of the time-dependent Schrodinger equation, 
which principle we had introduced almost three decades 
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ago yj4|. By this and with no further approximations or 
ad hoc assumptions, we derive a time-dependent effective 
potential with the following useful properties: (i) It satis- 
fies the exact large-separation asymptotic condition and 
is self-interaction free; (ii) It is expressed in terms of an 
equation easily solvable for both a finite and an infinite 
periodic (or homogeneous) problems; (iii) In the static 
case, our effective potential reduces to the earlier known 
localized Hartree-Fock (HF) potential [HI]. As an imme- 
diate application of this approach, we derive the dynamic 



exchange kernel f x {q, w) and the exchange shear modulus 
/i x of the homogeneous electron gas (HEG). 

Let us consider a A-electrons system with the Hamil- 
tonian 



m = E 



1a / 

--Aj + v ext (ri,t) 



E 



(i) 



The many-body wave-function ^(t) satisfies the 
Schrodinger equation 

i*(t) = H(t)*(t), 

where the dot denotes the time derivative. We are asking 
the question: What is the potential u e //(r, t) such that 
the functional 



t* a (t)-ir(t)* s (t) 



dri...dr 



N 



(2) 



is minimal at an arbitrary time t for the wave-function 
^ s {t) being the Slater determinant built with the single- 
particle orbitals tp a (r, t) which satisfy the single-particle 
Schrodinger equation 



iip a (r,t) = 



--A 



e//(r,t) 



(3) 



Let us assume that the TD part of v ext (r,t) was absent 
at t < 0, and we have already solved the static problem 
of determining v e f f (r) at t < and the corresponding 
orbitals ip a (r,t) &t t < 0. At t = 0, the time-dependence 
of the potential is switched on. Knowing the orbitals 
ip a (r,0), we find v e ff(r, 0) which, determining ip a (r,0) 
by Eq. ([3]), minimizes the functional ^ at t = 0. We 
then find ip a ( r i At) w "0a(r, 0) + V>a(r, 0)At, where At is 
a small time increase. The procedure is repeated up to 
an arbitrary time t. 

In the At — !• limit, this scheme reads: With fixed (but 
yet unknown) orbitals ip a (r, t), we are looking for the po- 
tential v e ff(r,t) which, determining ip a (r,t) by Eq. ©, 
minimizes the functional ([2]). This gives w e //(r,t) as a 
functional of the orbitals, and, finally, the orbitals them- 
selves are found by the self-consistent solution of Eqs. ((3|) . 
We note that a procedure of the minimization of the same 
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functional ([2]) with respect to ip a as independently varied 
functions retrieves the TD HF equations [l 21 ] - 

Since ^ s (t) obviously satisfies the equation 



4* 

2 



■V e ff(Ti,t) 



*.(*), 



(4) 



the functional ((2]) can be rewritten as 



l<7 



|* s (*)rdri...dr JV , (5) 



where S = v e ff — v ext . Equating to zero the first variation 
of Eq. <j5j) with respect to 5v = Sv, we find 



(6) 



<H> (r j , f ) dr i . . .dr jy = 0, 



which can be rewritten using the permutational symme- 
try of the wave-function as 



/ S>..*>-Ei^ 



|tf s (t)r<re(ri,i)dr 1 ...dr iv = 0, 



and. due to the arbitrariness of 8v, 



E fi ("'*)-Ei=rr 



ij - ij 



|*,(t)rdr 3 ...<frj V = 0. 



(7) 

Straightforward but rather lengthy transformations car- 
ried out in Ref. 15j lead from Eq. ([7]) to the following 
equation for the exchange potential v x = v — vh, where 
^^■(r) is the Hartree potential, 



n(r,t)v x (r,t) 



v x (ri,t)- 



1 



|p(r,n,t)| 2 dri 



r — ri 

p(r,r 1 ,t)p(r 1 ,r 2 ,t)p(r 2 ,r,t) 
|r x -r 2 | 



drirfr 2 , 



(8) 



where 



N 



n(r,t) = £|^(r,t)| 2 , (9) 

a=l 

N 

p(r,n,t) = £Va(r,*)C(ri,t) (10) 

are the particle density and the single-particle density- 
matrix, respectively [la ]. Equation (J5J is our main re- 
sult. We note that the only difference of Eq. ([8]) from 
the earlier known equation for the static localized HF 
potential (see Refs. 13 and [13 for spin- neutral and spin- 
polarized cases, respectively) is the time-dependence of 
all the quantities involved. It must, however, be em- 
phasized that without the derivation from the time- 
dependent variational principle, the generalization of the 
static localized HF potential to the time-dependent case 
by just inserting the time variable into the static equa- 
tion would have been ungrounded. Trivially, Eq. ([HJ re- 
duces to the equation for the localized HF potential in 
the time-independent case. 

Solution of Eq. flSJ in the case of a few-body system 
does not present a difficult problem, as have already been 
pointed out in Ref. [l3| in conjunction with the time- 
independent case: Due to Eq. (fTU|) . the kernel |p(r, ri, t)\ 2 
of the integral equation ((8]) is a separable function with 
respect to r and r' variables. Neither presents it a prob- 
lem in the case of infinite periodic systems, when the 
equation reduces to the matrix one. We now apply the 
TD localized HF potential to obtain the wave-vector and 
frequency dependent exchange kernel f x (q,uj) of HEG, 
which, being a fundamental quantity by itself, is also an 
important input in the theory of optical response of a 
weakly ingomogeneous interacting electron gas [18j . 

Dynamic exchange kernel of HEG - In the case of HEG 
and a weak externally applied potential Sv ex t(r,t) — 
3v ex t(c[,uj)e' l (' 1 ' r ~ ult \ we linearize Eq. (jHJ with respect to 
the latter potential. The zero-order orbitals are plane- 
waves [191 ] , and to the zeroth and first orders we have for 
the density-matrix 

p (r,r 1 ) = ^J2 I /( £ k)e 4k(r - ri) , 

k J 

6p(r, r 1 ,u) = Sv s (c L , u>) e *' *y JMzI&+sL e *<™i) 

V LO - 6k+q + £k + IT] 

(11) 



where ej. = fc 2 /2 are free-particle eigenenergies, /(et) are 
their occupation numbers, Sv s (q,uj) is the perturbation 
of the Kohn-Sham (KS) potential, V is the normaliza- 
tion volume, and r\ is an infinitesimal positive. After 
linearization, Eq. ([5]) yields 



3 



8v x (q,oj) _ Air 
6v s (q,u) {2irf 



1 



/(e k )B(g,k) - /(e k+q )[/(e k ) - l]C(fc) 
n - Ml) 



dk, (12) 



where 



S(q,k) 



(27T) 



2tt 



/(e k )/(e k+q )dk, 
-k x |2 



(2?r) 3 y |k 

2 /• /(e kl ) 



dki, 



k-kJ 



found as 

Sv x (q,uj) Sv s (q,uj) _ Sv x (q,uj) thy 1 / w j 
5« s (g,w) 6n(q,u) 6v s (q,u) 

where Xs (?) w ) i s Lindhard density- response function [2C 



2 ^ /(e k )-/(e k+q ) 



V ^ w-e k+q +e k + i?/ 



(17) 



In Ref. [13, we evaluate integrals (fT3"|) and (fl5)) analyti- 
cally and reduce the integral (ITU) to a single-fold one. 

It can be seen from Eqs. (|T2|) and (fTB"]) that Im/^ is 
nonzero inside the single particle-hole excitation contin- 
uum only. This deficiency of f x derived from TD local- 
ized HF potential is not limited to HEG, but, as can be 
easily seen, persists for any extended (periodic) system. 
Therefore, such subtle effect as the high-frequency tail 
of Jmf x of HEG dH cannot be accounted for within the 
present approach. Instead, we will now show that f x de- 
rived from the TD localized HF potential significantly 
corrects the Lindhard dielectric function of HEG within 
the single particle-hole continuum. 

In Fig. [TJ the exchange kernel f x (q,uj) obtained by the 
use of Eqs. (fT2"|) and (fTB|) is plotted at q = 0.5 x fc^ for 
r s = 2 and 5. With the inclusion of f x , the dielectric 
function of HEG can be written as 



-in 



q 2 1 + X?(«.W)/*(?,W) 



(18) 



In Fig. [2] the dielectric function of HEG of r s = 5 
obtained through Eq. (TT8|) is plotted together with the 
Lindhard dielectric function. From this we judge that 
dynamic exchange plays a significant role at this density 
and can hardly be considered as a weak perturbation. 

In Fig. 02 the static exchange kernel f x (q) is plotted 
for HEG of r s = 5 as a function of the wave-vector. Wc 
find a qualitative a gree ment with the Monte Carlo (MC) 
simulations of Ref. 22 . 
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FIG. 1. (color online) Exchange kernel f£(q, to), q = 0.5 x k,F, 
of HEG of r s = 2 (black curves online) and r a = 5 (red curves 
online). Solid and dashed curves are real and imaginary parts 
of/£. 




FIG. 2. (color online) Dielectric function e (q,u), q = 0.5 x 
kp, of HEG of r s — 5 evaluated by the use of EQ. (Tl8)) with 
the exchange kernel f£ included (red curves online) and its 
Lindhard counterpart el (q, oS) (black curves online). Solid and 
dashed curves are real and imaginary parts of the dielectric 
function, respectively. 



The following limiting cases can be further worked out 
from Eqs. (p ]) -(|T5 j) |23j 



(19) 
(20) 
(21) 
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4 




FIG. 3. (color online) Static exchange kernel f£ (q) of HEG of 
r a — 5. Solid line (red online) is the present result. Symbols 
with error bars are fx C (q) by MC simulations from Ref. [2^ |. 
The dashed line (blue online) shows the asymptotic behavior 
at large q, as stipulated by Eq. (|21[) . 



and since [14 1 



lim lim 

cj->0 g->0 



lim lim 

q— M uj — S-0 



4 ftxc 

3 n 2 



(22) 



where fj, x is the exchange contribution to the shear mod- 
ulus of HEG, we can write by virtue of Eqs. (|T^|) . (f2"U|) . 
and (HI 



3(1) 



3\2/3 



0.0091 



647r 5/3 r 4 



(23) 



which is smaller than the high-density result of Ref. [24 

3(f) 273 



/'rr 



0,0146 



407r 5/ 3r 4 — ■ In Table IH we compare \i 

obtained via Eq. §2?^ with [i xc of Refs. 0, [5 



and 24 



TABLE I. Exchange shear modulus ji x in units of 2u) p in. 
Present results are shown together with those for fi xc of Refs. 
GJ (QV),H (NCT), andlH (CV). 



present 0.01102 0.01559 0.01909 0.02204 0.02464 

QV 0.00738 0.00770 0.00801 0.00837 0.00851 

NCT 0.0064 0.052 0.0037 0.0020 0.0002 

CV 0.01763 0.02494 0.03054 0.03527 0.03943 

In conclusion, within the well-defined procedure of the 
minimization of the difference between the left- and right- 
hand sides of the time-dependent Shrodinger equation, 
we have derived a time-dependent single-particle effective 
potential for a system of arbitrary number of electrons 
under the action of a time-dependent external field. This 
potential is non-local, non-adiabatic, and self-interaction 
free, and it satisfies the exact large-separation asymp- 
totic condition. At the same time, our effective potential 



is comparatively easy for evaluation, which is in contrast 
to earlier known TD optimized effective potential. These 
properties open a way to efficiently use this potential 
within the context of time-dependent density-functional 
theory, as we demonstrate by the derivation of the ex- 
change kernel f x (q,uj) of the homogeneous electron gas. 

I acknowledge support from National Science Council, 
Taiwan, Grant No. 100-2112-M-001-025-MY3. 
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SUPPLEMENTAL MATERIAL 
TO THE PAPER BY V. U. NAZAROV 
"TIME-DEPENDENT LOCALIZED HARTREE-FOCK POTENTIAL" 

Derivation of Eq. fl8j. 

Using the permutational symmetry, we can rewrite Eq. ([7]) as 

(N-l) (N-l)(N-2) 



v(r u t) / |¥,(f)|"dr 2 ...dr 



N 



(N - l)v(r 2 ,t) - 



ri - r 2 



2|r 2 -r 3 | 



|* s (t)| 2 dr 2 ...dr w = 0, (S.l) 



z;(ri,t)n(ri,t) 

TV 



+ (iV-l) 



v(r 2 ,f) - 



N -2 



|ri-r 2 2|r 2 -r 3 | 



\V s (t)\ 2 dr 2 ...dr N =0, 



where n(r,t) is the particle-density. Further simplifications give 



7V(7V - 1)" 



ri - r 2 



n(ri,i) " 2 J n(r 1 , t)\r 2 - r 3 | 

where the fc-particle density-matrix is 

For the Slater determinant wave- function ^ s (t) the equalities hold 



p 2 (ri,r 2 ;ri,r 2 ;i) = 



N(N— 1) 
P3(ri,r 2 ,r 3 ;ri,r 2 ,r 3 ;i) = 



[n(ri,f)ra(r 2 ,i) - \p(r 1 ;r 2 ,t)\ 
1 



(S.2) 



ft.Cn, r 2 ; ri) r 2 ; t) ^ _ (iV-2) y ggC^£j^i^ dtadr8 = 0| (S . 3) 



(S.4) 



(S.5) 



7V(7V- l)(JV-2) 



[n(ri,t)n(r 2 ,t)n(r 3 ,t) - n(r 2 , t)|p(r x ; r 3 , i)| 2 - n(r 3 ,*)|/j(n;r 2 ,i)| 



- "-(i"i,i)|p(r 2 ;r 3 ,t)| 2 +p(ri;r 2 ,i) / 9(r 2 ;r 3 ,t)p(r 3 ;ri,t) +p(ri; r 3 , t)p(r 3 ; r 2 , i)/j(r 2 ; n, *)] , 



where 



p(r; r', t) = JV Pl (r; r', t) = £ ^(r, i)C (r', *)• 



(S.6) 
(S.7) 



Therefore 



v(r 2 ,t)-- 



where 



V H (r,t) 



n(r^) 

i / 1 ' 
r — r 



ri - r 2 



-/ 



n(ri,t) 7 n(ri,t)|r 2 - r 3 | 

/o(ri; r 2 , t)p(r 2 ; r 3 , i)p(r 3 ; n, t) 



n(ri,t)|r 2 - r 3 



dr 2 dr 3 + C(t) = 0, 



I -( *\ I +\a 1 f n(r 2 ,t)n(r 3 ,t) 1 f \p {r 2 ;r 3 ,t 2 

C(i) = / v(r 2 ,i)n(r 2 ,t)dr 2 - - / , : dr 2 dr 3 + - / — — dr 2 dr 3 = 0, 

2 J r 2 -r 3 2 J r 2 - r 3 



(S.8) 



(S.9) 
(S.10) 



and the second equality in Eq. (IS. 10[) follows from Eq. (|S.8[) . Eqs. (IS. 81) and (IS. 10[) give immediately Eq. ©. 



2 



Integrals (fl3)) -(fT5 j) . 



For integrals ([13]) and ([15]) we have straightforwardly 

1 



(2tt) 2 

where 0(x) is the Heaviside's step function and 



H(k,p) = - 



A{q) = ^©(2fcF - q)(2k F - g)'(4* F + g), 
= 77^H(k,k F ), 



(p 2 - k ) log 



p + k 



p — k 



2kp 



(S.ll) 
(S.12) 

(S.13) 



For integral ([14]) we have 



B(^) = j^®(k F - q )H(k,k F - q ) + J ^Q(2k F ~ q ) J P^k J 1 ^ q \ ^f)dk u 



(S.14) 



\k F -q\ 



where 



P(k,ki,x,y) = — { log 



+ 4fc 2 fc 2 a; 2 + 2fcfci (fcfci (2y 2 - 1) - 2a;y (fc 2 + fc 2 )) + kf + 2kk lX - y (fc 2 + fc 2 ) 



log((l-y)(fc-A ;i ) 2 )} 



(S.15) 



